Tracking solutions in tachyon cosmology 
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We perform a thorough phase-plane analysis of the flow defined by the equations of motion of a 
FRW universe filled with a tachyonic fluid plus a barotropic one. The tachyon potential is assumed 
to be of inverse square form, thus allowing for a two-dimensional autonomous system of equations. 
The Friedmann constraint, combined with a convenient choice of coordinates, renders the physical 
state compact. We find the fixed-point solutions, and discuss whether they represent attractors 
or not. The way the two fluids contribute at late-times to the fractional energy density depends 
on how fast the barotropic fluid redshifts. If it does it fast enough, the tachyonic fluid takes over 
at late times, but if the opposite happens, the situation will not be completely dominated by the 
barotropic fluid; instead there will be a residual non-negligible contribution from the tachyon subject 
to restrictions coming from nucleosynthesis. 



I. INTRODUCTION 

A phase of accelerated inflation in the early stages of 
our universe is favored by first-year WMAP data 0. 
Thus, it seems inflation is here to stay as the dominant 
paradigm for structure formation. The quest for a string 
theory motivated explanation of cosmological inflation 
has resulted in the emergence of the proposal of inflation 
driven by a tachyon field. The idea strongly relics in the 
possibility of describing tachyon condensates in terms of 
perfect fluids within string theories 0- A plethora of 
papers studying cosmological consequences of such fluids 
have appeared since, some in the framework of general 
relativity 0,0,1113, some others in the brane- world sce- 
nario 0- 

As happens with standard scalar fields, one's favorite 
inflationary behaviour is tailored by ad hoc choices of the 
initial conditions and the shape of the potential, but it is 
important to investigate up to what extent the features 
of the model depend on those choices. One way to ad- 
dress that problem is to consider tachyon field dynamics, 
because for a given potential such an analysis will provide 
us with constraints on the initial conditions. 

The stability of tachyonic inflation against changes 
in initial conditions has been studied for an exponen- 
tial potential and for the inverse power-law poten- 
tial |fj. Exact solutions for a purely tachyonic matter 
content with an inverse square potential are known 0, 
but no solutions exist for cases which combine tachyonic 
and barotropic fluids, so a dynamical systems approach 
may be relevant. Interestingly, the inverse square po- 
tential plays the same role for tachyon fields as the ex- 
ponenti al p otential |17| does for standard scalar fields 
U El EH E2] On the one hand, those are the po- 
tentials that give power-law solutions. On the other 
hand, only those potentials allow constructing a two- 
dimensional autonomous system |10| using the evolution 



equations, whereas for any other potential the number 
of dimensions will be higher if the system is to remain 
autonomous. 

As compared to [fj|, we throw in one more ingredi- 
ent in our study, because we allow for the presence of 
a barotropic fluid, along with the tachyon fluid. The 
crucial consequence is the appearance of fixed-point so- 
lutions in which the two fluids redshift at the same rate 
(tracking behaviour jjj), so that there is some sort of 
equilibrium. Tracking solutions are particularly inter- 
esting because their dynamical effects mimic a decaying 
cosmological constant (see 0, \^ for seminal refer- 
ences). Now, the fine-tuning problems posed by a cos- 
mological constant would be waived precisely because of 
the independence on the initial conditions. Nevertheless, 
the contribution of such relics to the fractional energy 
density are bounded by nucleosynthesis. 

On top of those interesting features, tracking solutions 
act as attractors at large, which means that the system 
is not any picky in what initial conditions are regarded. 

In Section II we study the phase-plane, find its fixed 
points and characterize them. In Section III we discuss 
the cosmological consequences of the attractor solutions: 
in subsection III A we consider tachyon dominated solu- 
tions, whereas in subsection III B we discuss the tracking 
ones. Finally, in Section IV wc outline our main conclu- 
sions and future prospects. 



II. PHASE PLANE 

The evolution equations for a flat (k = 0) Friedmann- 
Robertson- Walker (FRW) cosmological model filled with 
a tachyon field T evolving in a potential V(T) and a 
barotropic perfect fluid with equation of state p 7 = (7 — 
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Py + 3-fHpj = 0, (3) 
which arc, in turn, subject to the Friedmann constraint 
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Here and throughout overdots denote differentiation with 
respect to cosmic time t, H = a/a is the Hubble param- 
eter, and a is the synchronous scale factor. 

One can also define an energy density px and a pressure 
Pt for the tachyon, so that it can be thought of as a 
perfect fluid. We have then 
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If we use log a 3 as independent variable instead of the 
cosmological time, for any time dependent function / we 
get 



J 3H 

As usual, we also introduce n convenient variables: 
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Let us concentrate now on the inverse square potential 
V — j3T~ 2 . The evolution of the model is described by 
the dynamical system 



x' = (x 2 - l) (x - y/ay) , 

y' = V [x (x - y/ay) + z (7 - x 2 ) 

z' = z(z - 1) (7 - x 2 ) , 

along with the constraint 

, +2 = 1, 



(11) 

(12) 
(13) 



(14) 



which renders the phase space two dimensional so that 
we may speak of phase plane. For the sake of simplicity, 
we are using the following definition: 



(15) 



The physical constraints V > and p 7 > set limits 
on the dependent variables: —l<x<l,y>0 and 
x 2 +y 2 < 1. In consequence, the phase plane is restricted 
to the upper half of the unit disk centered at the origin, 
as depicted in Fig.^ For convenience we will also include 
the whole segment — l<x<l,y = 0in our phase plane. 




FIG. 1: Toroidal solenoid and point charge. 
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FIG. 2: Phase-space for a = 1.5 and 7 = 1.2. 

In addition, one can define a barotropic index for the 
tachyon fluid: 
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and, provided V ^ 0, one gets j T = T 2 . 

Depending on the values of 7 and j3 there may be 
up to five fixed points (O, A±, P and Q) and up to 
six heteroclinic orbits that connect pairs of fixed points 
(L±, C±, M±). The fixed point O located at the ori- 
gin (x, y) — (0, 0) corresponds to z = 1 and is a unsta- 
ble saddle (except in the very particular case in which 
7 = 0, which will be discussed below). The orbits L± 
in the stable manifold correspond to the characteristic 
exponent Ai = — 1 while the exponent for the unstable 
space is A2 = 7/2. 

The fractional densities of the two fluids are respec- 
tively defined as : 
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The fixed points A± located at (x,y) — (±1,0) corre- 
spond to z = 0, and are unstable nodes. The orbits L± 
in the unstable manifold correspond to the characteristic 
exponent Ai = 2 and orbits C± to A2 = 7/2. 

The fixed point P at (x, y) = (-Jctyi, yi), with 
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so that $la < for all < 7 < 71 , Q is always asymptot- 
ically stable and is a node (focus) when the argument of 
the square root is positive (negative). A particular case 
is shown in Fig. [31 



In the limit case in which 7 = 0, fixed points Q and O 
coincide and are the attractor in the system, as depicted 
in Fig. 2| 



FIG. 3: Phase-space for a — 1.5 and 7 = 0.3. 



always admits the arcs C± as the orbits in its stable 
manifold corresponding to the characteristic exponent 
Ai = — 1 + ay i/2 < 0, while the second exponent is 
A2 = ayi —7. In consequence, P is an asymptotically sta- 
ble node for 7 > 71 = ayi , in which case the phase space 
looks as depicted in Fig.El P is a global attractor, nearly 
all solutions ends there. When 7 = 71 a bifurcation 
arises: P turns into a unstable saddle and at the same 
time there appears a new attractor Q, which moves from 
P to O along an arc of the parabola (x, y) = (y/~f, j/ot) as 
7 decreases from 71 to (see Fig.^). The characteristic 
exponents are 
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FIG. 4: Phase-space for a = 1.5 and 7 = 0. 



TABLE I: The properties of the critical points. 

Name x y Existence Stability Ot 7t 

O All 7 and f3 Unstable saddle for 7 / Undefined 

Stable node for 7 = 

A+ 1 All 7 and (3 Unstable node 1 1 

A- - 1 All 7 and Unstable node 1 1 

P J/i All 7 and Stable node for 7 < 71 1 ayi 

Unstable saddle for 7 > 71 

7 7 

Q ^7 — 7 < 71 Stable node r^^= 7 

a fly 1 — 7 



III. COSMOLOGICAL FEATURES OF THE 
ATTRACTOR SOLUTIONS 



Now we turn our attention to the effective equation for 
the tachyon fluid in the attractor solutions, so that we 
can discuss in broad terms what sort of evolution they 
give rise to. 



A. Tachyon dominated solutions 

The attractor solutions (x,y) — (y^Syi, jvi) depict a 
situation in which the energy density of the fluid van- 
ishes, so it will be referred to as the tachyon dominated 
solution. It is straightforward to see it corresponds to 

7t = ayi (21) 
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and 

T = ^aylt + T Q , (22) 

with Tb an arbitrary integration constant. For the scale 
factor a we can set 

aoct 2/3Tr , (23) 

where the value of an integration constant has been fixed 
so that lim^o a = 0. 

The solution will be inflationary if px + ?>Pt < 0, and 
in terms of the tachyon field such condition is equivalent 
to 

T 2 < |, (24) 
which holds provided (3 > 2/ \fi. 

B. Tracking solutions 

We move on now to the most interesting case, which 
is represented by the attractor solutions with (x, y) — 
iVl,l/ a )- They satisfy 

Q T = 1 - 7 = 1 , (25) 

ay 1 - 7 

so that the energy density of the tachyon and of the 
fluid scale exactly as the same power of the scale factor, 
namely p 7 cx pt cx a -2 / 37 . For that reason, according 
to the definition by Liddle and Scherrer : 9|, these solu- 
tions display a tracking behaviour. Now, as in the case 
of the tachyon dominated solutions the expansion factor 
obeys a power-law, a oc i 2 / 37 , and inflation will proceed 
if 7 < min{2/3, ayi}. 

Nevertheless, apart from the considerations above, 
there exist restrictions on the values of a (or (3) that 
come from observations. Using standard nucleosynthe- 
sis and the observed abundances of primordial nuclides, 
the strong constraint that the fractional energy density 
of scalar matter cannot exceed 0.05 at temperatures near 
1 MeV was set [l4j . If we restrict the discussion to track- 
ing inflationary solutions (7 < 2/3), and set fir < 0.05, 



we see a > 23.09 is required. Such bound, though, can 
be evaded if the solution begins its history away from the 
tracking solution and in a region where fir -€i 1 (close to 
O) , and only reaches Q in the course of the evolution. 

IV. CONCLUSIONS 

The evolution equations of a spatially flat FRW uni- 
verse containing a barotropic fluid plus a tachyon T with 
an inverse-square potential V(T) = (3T~ 2 define a two- 
dimensional flow. The evolution of such models has been 
investigated by studying the orbits of that flow in the 
physical state, which is in this case a subset of the Eu- 
clidean plane. The Friedmann constraint, combined with 
a careful choice of coordinates, renders this subset com- 
pact. 

We have shown that the energy density of the tachyon 
dominates at late times for 7 > a(\/a 2 + 4 — a)/2, where 
7 is the barotropic index of the fluid and a = 4/3/3. In 
constraint, for 7 < a(V a 2 + 4 — a)/2, the barotropic 
fluid does not dominate completely and the contribution 
of tachyonic energy density to the total one is not neg- 
ligible. Nucleosynthesis imposes, then, tight bounds on 
the admissible values of a, but such restrictions can be 
relaxed if the locus of the initial solution is far from that 
of the tracking one and in a region where £It "C 1 (close 
to O), and only reaches Q in the course of the evolution. 

Finally, a possible generalization of this work would 
be considering generalized tachyon cosmologies like those 
presented in [lfij . 
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